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Abstract 

Irruptive forest insect pests cause considerable ecological and economic damage, 
and their outbreaks have been increasing in frequency and severity. We use a phy- 
logeographic approach to understand the location and progression of an outbreak 
by the MPB (Dendroctonus ponderosae Hopkins), an irruptive bark beetle that has 
caused unprecedented damage to lodgepole pine forests in western North America 
and is poised to expand its range across the boreal forest. We sampled MPB pop- 
ulations across British Columbia and Alberta and used phylogeographic methods 
to describe lineage diversification, characterize population structure, investigate ex- 
pansion dynamics, and identify source populations of the outbreak. Using 1181 bp 
of mitochondrial DNA sequence from 267 individuals, we found high haplotype 
diversity, low nucleotide diversity, and limited lineage diversification. The overall 
pattern was consistent with isolation by distance at a continental scale, and with 
reduced diversity and population structure in the northerly, outbreak regions. Post- 
Pleistocene expansion was detected, however more recent expansion signals were 
not detected, potentially due to the size and rapid rate of range expansion. Based 
on the limited genetic structure, there were likely multiple source populations in 
southern British Columbia, although the magnitude of the demographic expansion 
and rate of spread have obscured the signature of these source populations. Our data 
highlight the need for caution in interpreting phylogeographic results for species 
with similar demographics. 



Introduction 

The past few decades have seen an increase in the frequency 
and severity of insect pest outbreaks (Fleming and Candau 
1998; Bale etal. 2002; Jepsen etal. 2008), often impacting 
host species of economic importance. Understanding the dy- 
namics of an outbreak pest species is important from a man- 
agement perspective. Knowing the source and spread of an 
outbreak can help delineate potential boundaries for control 
and provide the basis for developing strategies to prevent 
future outbreaks (Mun et al. 2003; Kim et al. 2006; Porretta 
et al. 2007) . The use of population genetic methods to address 
management-associated questions for insect pest species is 
increasing. For example, Davies etal. (1999) determined that 
the source of a Californian outbreak of medfiy ( Ceratitis cap- 
itata) was Latin America. As well, Kobayashi etal. (2010) 
effectively used population genetic methods to determine 
that multiple source populations contributed to an outbreak 
of mirid bug (Stenotus rubrovittatus) in Japan, necessitating 



monitoring of nonoutbreak areas to prevent additional crop 
damage. Here, we use a phylogeographic approach to illumi- 
nate the outbreak dynamics of mountain pine beetle (MPB; 
Scolytidae: Dendroctonus ponderosae Hopkins, Fig. 1). 

The MPB is a bark beetle that attacks many Pinus species 
in western North America including: Pinus contorta Dougl, 
P. lambertiana Dougl., P. monticola Dougl. Ed. E. Don, P. 
ponderosae P. Laws. Ex C. Laws, P. albicaulis Engelm., and 
P. flexilis James (Furniss and Schenk 1969; Safranyik and 
Carroll 2006; Safranyik etal. 2010). Recently, MPB popula- 
tions have also been recorded in native stands of P. banksiana 
Lamb (Cullingham etal. 2011), an important boreal species 
in Canada and the eastern United States. While MPB are usu- 
ally considered secondary pests that feed on suppressed or 
diseased trees, certain conditions can lead to irruptive pop- 
ulation outbreaks. In these situations, MPB can successfully 
attack mature healthy trees, leading to significant stand mor- 
tality (Rudinsky 1962; Safranyik and Carroll 2006). These 
irruptive situations generally occur every 20-30 years and 
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Figure 1. Scanned election microscopy image of mountain pine beetle 
(Dendroctonus ponderosae). Photo credit: Jack Scott. 



usually last for — five to seven years (Cole and Cahill 1976; 
Safranyik 2004; Taylor etal. 2006). 

Outbreaks of MPB are important for the maintenance of 
healthy forest stands (Roe and Amman 1 970; Logan and Pow- 
ell 2001; Taylor and Carroll 2004), but anthropogenic factors 
have altered the dynamics of this system. Increasing global 
temperatures have prompted synchrony of beetle emergence 
in the summer and reduced winter mortality, thereby result- 
ing in population increases (Hicke et al. 2006; Regniere and 
Bentz 2007; Powell and Bentz 2009). Population growth of 
MPB has also been fueled by an abundant food supply of 
continuous tracts of even-aged forest stands resulting from 
forest management practices (Taylor et al. 2006). All of these 
factors have contributed to the most severe MPB outbreak on 
record, and an unprecedented eastern expansion into north- 
ern Alberta outside of the historical range (Taylor and Carroll 
2004; Raffa et al. 2008; Robertson et al. 2009; Bentz et al. 2010; 
Safranyik etal. 2010; Cullingham etal. 2011). 

As MPB populations expand into novel habitats, we have 
been presented with an opportunity to explore their phy- 
logeographic patterns in real time. Understanding how and 
where spatial expansion has occurred can assist with MPB 
management and help to predict risk of future spread. Phy- 
logenetic analysis has been previously employed to analyze 
MPB populations in western Canada to understand historic 
demographics (Mock etal. 2007), although sampling was 
diffuse across the entire range, limiting the inferences that 
could be drawn about MPB phylogeography. Increased sam- 
pling within the outbreak populations would provide a more 
comprehensive understanding of the outbreak history. We 
are provided with a unique opportunity to apply expan- 
sion statistics to an outbreaking natural system. The range 
of MPB was affected by the Pleistocene glaciation due to 
range contraction and expansion of its host species caused by 



glaciation (Fazekas and Yeh 2006) . In addition to this historic 
expansion, their range has grown considerably over the past 
50 years into northern BC and west-central Alberta, repre- 
senting a recent expansion. Here, we analyze MPB popu- 
lations in British Columbia and Alberta, focusing on these 
areas of range expansion. We have included mitochondrial 
(mtDNA) data from Mock et al. (2007) to assess range-wide 
lineage diversification. The objectives of our study were to ( 1 ) 
assess lineage diversification and patterns of diversity in re- 
lation to geography across the MPB range, (2) to determine 
if signals of expansion can be detected in a natural system 
with complex population dynamics comprising historic and 
recent expansions, and (3) determine whether these data can 
be used to identify locations of source populations that may 
have contributed to the current outbreak. 

Methods 
Sample collection 

Live beetles were collected from 26 locations in British 
Columbia, Alberta, and South Dakota (Table 1 and Fig. 2) 
from 2005 to 2009 prior to summer dispersal each year. Bee- 
tle larvae and adults were collected from different galleries by 
direct sampling from infested lodgepole pine and lodgepole 
pine x jack pine hybrid trees. For each tree, the Global Posi- 
tioning System location was taken and collected beetles were 
transported to the laboratory on ice or in 95% ethanol and 
stored at -80° C upon arrival. We have also included in our 
analyses MPB samples analyzed by Mock et al. (2007) at eight 
localities: Arizona (F-AZ), British Columbia (FSJ-BC), Cali- 
fornia (R-CA, SB-CA), Idaho (BF-ID, S-ID), Oregon (L-OR), 
and Utah (K-UT) (GenBank: DQ865977. 1-DQ86602 1 . 1 ; Ta- 
ble 1 and Fig. 2). 

DNA extraction and sequencing 

Genomic DNA was isolated from whole beetle samples using 
one of two methods. The first method employed a standard 
phenol/chloroform procedure (Sambrook and Russell 2001) 
and following precipitation, DNA was suspended in Tris- 
EDTA (pH 8.0). Alternatively, a DNeasy Blood & Tissue Kit 
(Qiagen, Mississauga, Ontario) was used following manufac- 
turer's instructions. 

Partial mtDNA sequence was generated from portions of 
cytochrome c oxidase I (COI) and cytochrome c oxidase II 
(COII) genes as well as the complete tRNA-LEU using two 
sets of primers (Cl-J-1718 (F): GGA GGA TTT GGA AAT 
TGA TTA GTT CC, Cl-N-2611 (R): GCA AAA ACT GCA 
CCT ATT GA and Cl-J-2193 (F): CCA GGG TTT GGT ATA 
ATT TCT T, C2-N-3174 (R): TTA GAG GGG AAG ACC TAT 
CTT GT; Simon et al. 1994; Cognato and Sperling 2000; Mock 
etal. 2007). PCR amplification was carried out in 25-ul re- 
actions containing 25 ng DNA, 1 x Thermopol buffer (NEB, 
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Table 1. Sample areas for Dendroctonus ponderosae mitochondrial DNA analysis with sample size (W), haplotype diversity (H), nucleotide diversity (n) 
(and associated standard deviation "SD"), their population assignment based population histories, and associated diversity measures. Sample areas in 
bold are taken from Mock etal. (2007). 



Pop 


Abb rev. 


Year 


W 


H 


SD 


71 


SD 


North 






173 


0.9035 


0.0109 


0.0065 


0.0034 


Bowron Lake 


BL 


2006 


13 


0.8590 


0.0886 


0.0087 


0.0048 


Fox Creek 


FC 


2008 


12 


0.9242 


0.0575 


0.0057 


0.0033 


Francois Lake 


FL 


2006 


1 1 


0.8909 


0.0740 


0.0029 


0.0018 


Fairview 


FV 


2008 


7 


0.9048 


0.1033 


0.0082 


0.0049 


Grande Prairie 


GP 


2008 


10 


0.8444 


0.0796 


0.0044 


0.0027 


Houston 


HO 


2006 


12 


0.8636 


0.0639 


0.0014 


0.0010 


McBride 


MB 


2005 


10 


0.9778 


0.0540 


0.0083 


0.0047 


Mackenzie 


MK 


2005 


9 


0.8889 


0.0910 


0.0089 


0.0051 


Mount Robson 


MR 


2005 


10 


0.9556 


0.0594 


0.0085 


0.0048 


Prince George 


PG 


2005 


9 


0.9444 


0.0702 


0.0068 


0.0040 


Pine Pass 


PP 


2006 


7 


0.8095 


0.1298 


0.0047 


0.0029 


Quesnel 


QN 


2006 


7 


0.9048 


0.1033 


0.0101 


0.0060 


Telkwa 


TK 


2006 


9 


0.8889 


0.0910 


0.0058 


0.0034 


Tumbler Ridge 


TR 


2007 


17 


0.8309 


0.0648 


0.0035 


0.0020 


Valemount 


V 


2007 


10 


0.8444 


0.1029 


0.0089 


0.0050 


Willmore-Kakwa 


WK 


2008 


11 


0.9818 


0.4630 


0.0091 


0.0051 


Ft St James 


FSJ-BC 


2003 


9 


0.2222 


0.1662 


0.0022 


0.0015 


Central 






78 


0.9624 


0.0106 


0.0089 


0.0045 


Cypress Hills 


CH 


2007 


12 


0.9545 


0.0569 


0.0115 


0.0063 


Canmore 


CA 


07/08 


16 


0.9250 


0.0389 


0.0104 


0.0056 


Crowsnest Pass 


CN 


2008 


12 


0.9545 


0.0569 


0.0052 


0.0030 


Golden 


G 


2008 


7 


0.9524 


0.0955 


0.0109 


0.0067 


Kotonay-Yoho 


KY 


07/08 


6 


1.0000 


0.0962 


0.0094 


0.0058 


Sparwood 


SW 


2008 


8 


0.9643 


0.0772 


0.0082 


0.0048 


Whistler 


WH 


2006 


8 


1.0000 


0.0625 


0.0070 


0.0042 


Bonner's Ferry 


BF-ID 


<2006 


9 


1.0000 


0.0524 


0.0694 


0.0375 


South 






87 


0 9588 


0 0096 


0 0087 


0 0045 


Black Hills 


BH 


2009 


11 


0.8545 


0.0852 


0.0085 


0.0047 


Custer Peak 


CU 


2009 


10 


0.7556 


0.1295 


0.0054 


0.0032 


Tigervill 


TV 


2009 


13 


0.8462 


0.0854 


0.0054 


0.0031 


Flagstaff 


F-AZ 


<2006 


10 


0.5333 


0.0947 


0.0108 


0.0060 


Kamas 


K-UT 


<2006 


10 


0.7778 


0.1374 


0.0849 


0.0452 


La Grande 


L-OR 


<2006 


8 


1.0000 


0.0625 


0.1082 


0.0594 


Klamath 


R-CA 


<2006 


10 


0.9556 


0.0594 


0.0834 


0.0444 


San Bernadino 


SB-CA 


<2006 


8 


0.7500 


0.1391 


0.0897 


0.0492 


Stanley 


SID 


<2006 


7 


0.7619 


0.1148 


0.0846 


0.0476 



Pickering, Ontario), 0.1 ng/jul bovine serum albumin frac- 
tion V (BSA), 4 mM MgCl 2 , 0.2 mM each dNTP, 0.25 /iM 
each primer, and 1U Taq DNA polymerase (NEB). Thermo- 
cycling conditions consisted of an initial denaturation at 98°C 
for 30 sec, followed by 35 cycles of 95°C for 10 sec, 45°C for 
30 sec, and 72°C for 60 sec, and a final extension at 72.0°C for 
5 min. Forprimer pair Cl-J-2193 and C2-N-3174, 15 ngDNA 
was used and the annealing temperature was 62° C. Primers 
and unincorporated dNTPs were removed using ExoSAP-IT 
(USB Corporation, Cleveland, OH) following manufacturer's 
instructions prior to sequencing. Purified PCR products were 
sequenced bidirectionally with Applied Biosystems Big Dye 
Terminator v 3.1 cycle sequencing kit (Applied Biosystems, 



Foster City, CA) using an Applied Biosystems 3730 DNA 
Analyzer. Sequence data were analyzed and assembled with 
Seqman (DNASTAR Lasergene, Madison, WI). Final assem- 
bled sequences were checked manually. 

Phylogenetic analyses 

The most appropriate model of sequence evolution was 
selected using jModelTest 0.1.1 (Posada 2008) based on 
Akaike's information criterion. Using the parameter esti- 
mates, we calculated both a Maximum-Likelihood (ML) 
tree and a Median-Joining (MJ) network for the aligned 
mtDNA dataset. The ML tree was estimated using the on-line 
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Figure 2. Sample sites, including data from Mock etal. (2007) divided into south (blue), central (green), and northern (red) groups. Beetle attack data 
for British Columbia and Alberta up to 1970 are included as shading as an indicator of the historic range postglaciation in Canada. The haplotype 
Median-Joining network for D. ponderosae is included, where the presence of a haplotype is indicated for each group. Circles are proportioned based 
on sample size. Sample sites in Arizona and California from Mock etal. (2007) are not shown. Beetle attack data are from Thandi and Taylor (unpubl. 
ms.). 



version of PhyML 3.0 (http://www.atgc-montpellier.fr/ 
phyml/; Guindon and Gascuel 2003) using a general time 
reversible (GTR) model of nucleotide substitution, with the 
transitiomtransversion ratio, invariable sites, and gamma- 
shape parameter estimated from jModelTest, subtree prun- 
ing and regrafting for tree improvement with 10 starting 
trees, Shimodaira-Hasegawa like branch support, and 100 
bootstrap replicates. We included D. jeffreyi Hopkins, sister 
species to D. ponderosae (Kelley and Farrell 1998) as an out- 
group. The MJ network was calculated using Network 4.516 
(fluxus-engineering.com). We implemented the maximum 



parsimony option to identify unnecessary median vectors 
and links (Polzin and Daneschmand 2003) and calculated 
the MJ network (Bandelt etal. 1999) using the output from 
this and the parameters from jModelTest. 

Using a traditional phylogeographic approach (Avise 
2004), we looked at lineage diversification in relation to 
geography across the range of the beetle. We divided our 
study area into three separate regions based on their bio- 
geographic history (Fig. 2): (1) south — populations that 
were south of the Cordilleran ice sheet (Clark etal. 2009), 
(2) central — populations that are considered to have been 
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endemic since post-Pleistocene expansions (Safranyik etal. 
2010), (3) north — populations that have resulted from a re- 
cent range expansion of MPB (Fig. 2). 

Population structure 

Given the evidence of post-Pleistocene colonization (Mock 
et al. 2007) and a more recent expansion, we calculated hap- 
lotype (H) and nucleotide (jr) diversity for each subpopu- 
lation, and regressed these on latitude and longitude, where 
we expected diversity to decrease toward the north and to 
the east as MPB has recently expanded into central Alberta. 
All above calculations were performed in Arlequin 3.5.1.2 
(Excoffier etal. 2005) and significance was estimated using 
10,000 bootstrap replicates where appropriate. To test for sig- 
nals of expansion, we estimated Fu's Fs (Fu 1997), and R2 
from Ramos-Onsins and Rozas (2002) given the statistical 



power of these statistics over other measures of expansion 
(Ramos-Onsins and Rozas 2002). Both measures were calcu- 
lated in DnaSP v 5 (Librado and Rozas 2009) and significance 
was tested using the coalescent simulation to estimate the 
distribution. Using continuous surface maps, we visualized 
geographic patterns of diversity and signatures of expansion. 
The maps were generated in ArcMap 9.2 (ESRI) using the 
inverse-distance weighting (Fig. 3; Shepard 1968) calculation 
in the spatial analyst extension, using the six nearest neigh- 
bors for the estimation. 

To examine the spatial genetic structure of MPB pop- 
ulations, we tested for isolation by distance including the 
Mock etal. (2007) data, and for just the outbreaking re- 
gion in British Columbia and Alberta. We calculated pair- 
wise genetic differentiation using F ST in Arlequin, which ac- 
counts only for haplotype frequency differences (Excoffier 
etal. 1992). Geographic distances were calculated among 
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sample area centroids using Hawth's Tools (Beyer 2004) in 
ArcMap 9.2 (ESRI Inc, Redlands, CA). The correlation be- 
tween the matrices was tested using a Mantel test imple- 
mented in zt (Bonnet and Van de Peer 2002) and significance 
was assessed using 10,000 iterations. 

Results 

Phylogenetic analyses 

We obtained 1 18 1 bp of mtDNA sequence from 267 individu- 
als yielding 62 haplotypes (GenBank: JQ308436-JQ308479). 
We aligned our sequences with those from Mock et al. (2007) 
trimmed to our fragment size. This reduced the number of 
haplotypes described in Mock et al. (2007) to 35, which were 
analyzed with the 62 haplotypes from the newly sampled 
MPB populations. Most of the haplotypes were found in one 
of the three regions (south [31], central [23], andnorth [26]), 
1 1 were shared between two of the regions (south-central 
[6], south-north [3], and central-north [2]), and six were 
shared among all regions. One hundred six variable sites were 
found, 97 transitions, seven transversions, and two transi- 
tion/transversions. The ML tree was very shallow, character- 
ized by short internal branches with little bootstrap support 
over 80% (Fig. SI). The low level of divergence among hap- 
lotypes was also reflected in the MJ network where there 
were some reticulations and a number of "star-burst" pat- 
terns (Fig. 2). Some phylogenetic resolution in both the MJ 
network and ML tree was found among haplotypes that oc- 
curred in the southern portion of the MBP range, although 
for the majority, lineage diversification did not correspond 
to geography (Fig. 2). 

Population structure 

Among sampling areas, we found high haplotype diversity 
(0.756-1.000) and low nucleotide diversity (0.001-0.011) 
(Table 1). Nucleotide diversity declined significantly with in- 
creasing latitude (R 2 = 0.39, t 32 = -3.41 x 10~ 3 , P = 1 Al x 
1 0 -5 ) , but haplotype diversity was unrelated to either latitude 
or longitude (£ 2 i at itude = 0.05, t 32 = 5.94 x 10 3 , P = 0.19; 
^longitude = 1-38 x 10" 3 , f 32 = -9.42 x 10~ 5 , P = 0.83). 
There was a greater than 10-fold difference in nucleotide di- 
versity between the southern and northern beetle populations 
(Fig. 2). There were minor differences between the two ex- 
pansion statistics we calculated; therefore, we chose to map fs 
to illustrate patterns of demographic expansion (Fig. 4) . Only 
two populations (BF-ID and WH) had significant signals of 
expansion, located in the central regions, we plotted the ac- 
tual values of Fs because they were highly correlated with the 
P- values (r = 0.93). The signal of expansion decreased to the 
north with no evidence of expansion in the south. However, 
when we analyzed the population groups based on their his- 
tories, we found a significant signal of expansion only for the 



central (Fs = -13.40, P = 0.003) and no significant signal 
for the southern (Fs = -7.79, P = 0.050) or northern groups 
(f s = -5.76, P = 0.107). Genetic differentiation was highly 
correlated with geographic distance for the combined data 
(r = 0.64, P = 0.0001, Fig. 5A), however this relationship 
was weaker when only the current outbreak samples were 
included (r = 0.30, P = 0.0033, Fig. 5B). There were no sig- 
nificant F ST comparisons among our 22 sampled populations 
in the outbreak region following Bonferroni correction. 

Discussion 

Our assessment of mitochondrial diversity in MPB across 
British Columbia and Alberta indicates poorly resolved 
phylogenetic structure with limited genetic differentiation. 
Across the sampled range, we observed high haplotype diver- 
sity, but found over a 10-fold difference in nucleotide diversity 
between the southern and central/northern populations, re- 
sulting in a significant correlation with latitude. When com- 
bined with the data from Mock etal. (2007), we found a 
significant signal of isolation by distance, which agrees with 
their prior analysis. Based on the limited genetic structure, we 
hypothesize there were multiple outbreak sources in south- 
ern British Columbia in the current epidemic. Analysis with 
microsatellite markers also found southern British Columbia 
represented a single genetic population (N. G. Samarasekara 
etal. unpubl. ms.) with weak substructure due to the size 
and rapid spread of this outbreak. Although little phyloge- 
netic structure was apparent, there were interesting patterns 
of diversity and expansion that can be used to make some 
demographic inferences. As well, these data bring to light po- 
tential limitations of phylogeographic inference for species 
with similar demographics. 

Limited phylogenetic structure of MPB is not surprising 
given their population demographics. MPB development is 
completed most often in one year and under certain condi- 
tions in two years (Amman 1977; Safranyik etal. 2010) and 
following development within the host tree, adult beetles 
emerge and disperse (Safranyik and Carroll 2006). Disper- 
sal is either short (among host trees, 30 m; Safranyik etal. 
1992) or medium distance (between forest stands, 20 km; 
Furniss and Furniss 1972), although beetles have been shown 
to achieve long-range dispersal of >100 km in a single day 
(Safranyik etal. 2010). Given the dispersive nature of the 
MPB life cycle and the distances they travel, the potential for 
isolated populations and genetic signatures to develop should 
be limited. 

Along with being dispersive, MPB also routinely experi- 
ence large fluctuations in population size. Since documen- 
tation of MPB attacks began in British Columbia approxi- 
mately 125 years ago, there have been four to five major out- 
breaks, often in different geographic areas (Taylor and Carroll 
2004; Aukema etal. 2006). The level of haplotype diversity 
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Figure 4. Spatially interpolated surface of range expansion (measured as Fu's [1 997] Fs) for MPB populations across western North America, generated 
using inverse-distance weighting in the spatial analyst extension of ArcMap9.2. 



is affected by the effective population size, where a loss of 
diversity will result when the effective population is reduced. 
For outbreaking populations, there is an expectation of re- 
duced diversity since the long-term effective population will 
reflect the size of the endemic phase (Motro and Thomson 
1982). Estimates of population size for MPB in the endemic 
phase are <40 attacking females per hectare (Carroll etal. 
2006), a considerably smaller population than the epidemic 
phase where attack densities are on average 70 beetles/m 2 
bark (Raffa and Berryman 1983). Despite this expectation, 
we observed very high haplotype diversity in all of our studied 
populations (Table 1), consistent with other species of irrup- 
tive bark beetles (Cognato etal. 1999; Anducho-Reyes etal. 



2008; Ruiz et al. 2009). Based on other studies of outbreaking 
populations, there are two factors that likely act to maintain 
high diversity. First, the dispersive behavior of MPB will re- 
sult in constant immigration, which can act to maintain high 
levels of heterozygosity (Ehrich and Jorde 2005). Second, the 
effective population size at the endemic phase for MPB may 
be high enough to maintain diversity, and/or the length of 
the endemic phase is not long enough for a considerable de- 
crease in diversity, this was observed for outbreaking locust 
populations in Europe (Chapuis etal. 2009). 

Using a phylogeographic approach, other studies have been 
able to assess the geographic source of outbreaking insects 
(Davies etal. 1999; Kobayashi etal. 2010). When multiple 
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Figure 5. Isolation by distance plot for D. ponderosae across (A) western North America; and (B) our sample areas in British Columbia and Alberta. 
Geographic distance is significantly correlated to genetic differentiation (F ST ) across the entire region (r = 0.64, P = 0.0001), though more weakly 
among the outbreak populations alone (r = 0.30, P — 0.0033). Comparisons with FSJ-BC are indicated by gray squares. 



sources contribute to an outbreak, it will generally result in 
genetic population clusters (Mun et al. 2003; Kobayashi et al. 
2010). Spatiotemporal analysis of outbreak data in BC in- 
dicates the current epidemic resulted from a simultaneous 
increase in populations in southern BC (Aukema et al. 2006), 
suggesting multiple sources. Yet, we did not observe evidence 
of population structure, rather we observed nonsignificant 
pair- wise F S t comparisons and very diffuse isolation by dis- 
tance among outbreak populations (Fig. 5B). In both Mun 
etal. (2003) and Kobayashi etal. (2010), the outbreaks were 
geographically separate whereas the MPB outbreak in south- 
ern BC involved multiple areas that are geographically prox- 
imate. Proximity and the dispersal capability of MPB have 



likely contributed to the mixing of mitochondrial lineages 
during the outbreak, which would be consistent with the 
findings of Aukema etal. (2006). 

The limited degree of phylogenetic structure across a large 
geographic region for MPB raises an important considera- 
tion: how applicable is a phylogeographic approach to re- 
solving the history of a highly dispersive, irruptive species? 
There have been a number of studies that have considered the 
phylogeographic structure of outbreaking insect populations 
(Cognato etal. 2005; Maroja etal. 2007; Mock etal. 2007; 
Anducho-Reyes et al. 2008; Ruiz et al. 2009) and these found 
limited phylogenetic structure, high haplotype diversity, low 
nucleotide diversity, and evidence of range expansion. For 
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both our analysis of MPB and that of Mock etal. (2007), the 
signal of range expansion that we observe is most likely to 
be from post-Pleistocene historical processes and the more 
recent expansion is not detectable using traditional expan- 
sion statistics. Since record taking began in the early 1900s, 
MPB has been documented primarily in southern and cen- 
tral British Columbia (Carroll etal. 2003; Taylor and Carroll 
2004), and in the past 50 years the range has increased sub- 
stantially toward the north. While we were unable to detect 
expansion using statistical measures, we do see the relation- 
ship between haplotype diversity and latitude that is sugges- 
tive of a historical range expansion. 

The magnitude of the recent outbreak, in terms of both the 
rapid population size increase and the geographic spread, has 
the potential to swamp phylogeographic signal. For example, 
Mock etal. (2007) sampled one site in British Columbia at 
Fort St. James (FSJ) in 2003 (S. Lindgren, pers. comm.). FSJ 
is proximate to our sample sites in Francois Lake (FL) and 
Mackenzie (MK), which were obtained in 2005/2006. The FSJ 
sample is characterized by extremely low haplotype diversity 
(H = 0.2222) relative to all of our datapoints including its 
nearest neighbors (H FL = 0.8909, H MK = 0.8889). This low 
diversity is a signature of a very recently founded population 
(Mayr 1942), and in 2003 this area was at the front of the 
MPB outbreak. Yet, within three years (equal to three genera- 
tions), this signature is no longer evident in our samples. This 
suggests that recent population demographic change has the 
potential to obscure population genetic signatures. In fact, 
Lee et al. (2007) studied the genetic diversity of Grosmania 
clavigera, a fungal associate of MPB, and found two distinct 
lineages in the northern Rocky Mountains. However, more 
recent studies did not find any lineage diversification in a 
comprehensive treatment of the region (Alamo uti et al. 201 1; 
Roe etal. 2011), suggesting that rapid expansion of MPB 
replaced the endemic populations of fungi in the northern 
Rocky Mountains. 

This study provides information on MPB dynamics, but 
we have also highlighted potential drawbacks in using phy- 
logeographic analyses for highly dispersive species with large 
fluctuations in population size. This is clearly an avenue that 
deserves more detailed attention through model-based phy- 
logeographic inference (Excoffier and Heckel 2006; Kuhner 
2008; Beaumont et al. 2010) as the use of phylogeography in 
these types of applications increases (Mun et al. 2003; Cog- 
nate et al. 2005; Kim et al. 2006; Maroja et al. 2007; Mocket al. 
2007; Porretta etal. 2007; Anducho-Reyes etal. 2008; Cha- 
puis etal. 2009; Ruiz etal. 2009). The use of a single marker 
may affect our ability to detect accurate patterns, a common 
criticism of single marker studies (e.g., Zink and Barrow- 
clough 2008; Flanders etal. 2009); however, multiple nuclear 
markers are not likely to contradict these results, particularly 
in outbreak areas that are phylogeographically unstructured. 
A recent population genetic analysis using microsatellites did 



not reveal conflicting or additional information about MPB 
population history (N. G. Samarasakara etal. unpubl. ms.), 
rather the higher resolution of microsatellite data was able to 
resolve subtle population structure resulting from the recent 
expansion (N. G. Samarasakera etal. unpubl. ms.). Also, mi- 
crosatellite data may not be appropriate for phylogeographic 
analysis given the differences in temporal sensitivity of mi- 
crosatellite and mitochondrial data (Wang 2010). 
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